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Abstract 

Probing the ss content of the n and n' mesons and considering mixing between these states 
as well as gluonic contributions, the form factors responsible for semileptonic D s — > (ry, n')lv 
transitions are calculated via light cone QCD sum rules. Corresponding branching fractions and 
their ratio for different mixing angles are also obtained. Our results are in a good consistency 
with experimental data as well as predictions of other nonperturbative approaches. 
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I. INTRODUCTION 



Based on experimental results, a considerable part of the total decay rate of the D s 
meson is related to its decay to rj and rj' mesons. Therefore, the D s is a proper meson to 
study the phenomenology of the r\ and r/ mesons and their structures. Due to charm quark, 
this meson plays an essential role in analyzing of the weak and strong interactions as well as 
exploring new physics beyond the standard model (SM) which will be probed by the large 
hadron collider (LHC). The charmed systems are known for very small CP violations in the 
SM, hence any detection of CP violations in such systems can be considered as a signal for 
presence of new physics (for more information about the D s meson and its decays see 

In the present work, we analyze the semileptonic D s — > (77, rf)lv decays in the framework 



of light cone QCD sum rules (LCSR). The rj and rf mesons are mixing states 



A3, 



\r]) = cos ip\r) q ) - sin <p\r) s ), 

\rf) = sxa cp\r] q } + cos (p\r] s ), (1) 

where tp is single mixing angle. The measured values of ip in the the quark flavor (QF) 
basis (for more information about this basis see for instance ^-7]) are (p = (39.7 ±0.7)° and 
41.5 d= 0.3 sta t ± 0.7 syst ± 0.6 t h)° with and without the gluonium content for rj', respectively 
sj]. The mixing angle ip has also been obtained as <p — [39.9 ± 2.6(exp) ± 2.3(th)]° by 
recently measured BR[D(D S ) — > 77(7/) + I + v{\ in light-front quark model 
In QF basis, 

\Vq) — —7= (\uu) + \dd) 



V2 

\Vs) = \ss) . (2) 

Since the D s meson decays to rj and 7/ via rj s state, the transition form factors of these 
decays in the QF basis are written in terms of the transition form factors of D s — > r/ s as: 

= ~ sm <p x , f^' = cos <p x . (3) 

For calculation of /, Ds_i-r ' < ) via the LCSR through f® 3 ^ 3 , information about distribution 
amplitudes (DA's) of the \rj s ) state as well as corresponding parameters are needed. These 
quantities have not been known yet, exactly. However, the same quantities for rj meson are 



2 



available and investigation of f^ s ~* v is possible, directly. On the other hand according to 
Eq. (EJ), there is a relation between f^ s ^ n and f[ >s ~' ri , 

— n — , = tano9, 4 

so, our strategy will be as follow. First, we will calculate the form factors, f® 3 ^ via 
the LCSR, then using Eq. (j4]) and the values of the mixing angle ip, we will evaluate the 
transition form factors of D s — > r(lv. 

The paper is organized as follows. In the next section, we obtain the LCSR for the 
transition form factors responsible for D s — > rjlu decay. Section III is devoted to the 
numerical analysis of the form factors and calculation of branching ratios of the D s — > 
(rj,ri')lv decays. We also compare the obtained results with the existing predictions of the 
other nonperturbative approaches as well as experimental data. 



II. LCSR FOR D s -> 7] TRANSITION FORM FACTORS 

To calculate the transition form factors of the D s — > r] in LCSR method, we consider 
the following correlation function: 

n fM (p,q)=ij d 4 a;e^(r / (p)|T{^)7 M (l-75)c(x)c(0)z(l- 75)5(0)} |0), (5) 

where we will use the DA's of the t] meson. The main reason for choosing the Chiral 
current, ci{l — 7 5 )s instead of the usual pseudoscalar (PS), ci^s is to eliminate effectively 
the contribution of the twist-3 wave functions which are poorly known and cause the main 
uncertainties to the sum rules. This current provides results with less uncertainties (see 



also 



10 



14|). Here, we should stress that the Chiral current may enhance the NLO twist-2 
contribution and to get more exact results, one should use the DA's of the r\ mesons up to 
NLO which are not available yet. 

According to the general philosophy of the QCD sum rules and its extension, light cone 
sum rules, we should calculate the above correlation function in two different ways. In phe- 
nomenological or physical representation, it is calculated in terms of hadronic parameters. 
In QCD side, it is obtained in terms of DA's and QCD degrees of freedom. LCSR sum 
rules for the physical quantities like form factors are acquired equating coefficient of the 



3 



sufficient structures from both representations of the same correlation function through dis- 
persion relation and applying Borel transformation and continuum subtraction to suppress 
the contributions of the higher states and continuum. 

To obtain the phenomenological representation of the correlation function, we insert a 
complete set of D s states between the currents. Isolating the pole term of the lowest PS 
D s meson, we get, 



n (« g ) = (^(P)l^(l - l5)c\D s (p + q))(Ds(p + gMl - 7 5 )s|0) + _ (g) 



where ■ • • stands for contributions of the higher states and continuum. The matrix element, 
(D s \ci(l — 75)s|0) is defined as: 

( J D s |cz(l- 75 ) S |0) = ^%^, (7) 

where fjj a is leptonic decay constant of D s meson. The transition matrix element, 



siderations as 
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J2s(p + q)) can be parameterized via Lorentz invariance and parity con- 



12: 



(i7(p)|5T M (l - ls)c\D s (p + q)) = 2f^(q 2 )p, + (f^(q 2 ) + f^\q 2 ))q, , (8) 

where, f± B ^ v (q 2 ) are transition form factors responsible for D s — > r) decay. Using Eqs. (j7j) 
and (jSj) in Eq. we obtain, 

n„(p, q) = Il^q 2 , (p + q) 2 )p, + Il 2 (g 2 , (p + qf)q^ (9) 

where, 

2f^(q 2 )m 2 f Ds r ^ (s) 



ITi = t T7 — 5 7 — s N0N + / ds — + subtractions , 

{m c + m s )(m 2 Dg - {p + q) 2 ) Js s - (p + q) 2 

n (tf'^(g 2 ) + /^V)K>>. r j pM u + t - 

II2 = — 7 T7 — 5 7 vfr 1- / ds -. + subtractions , (10) 

[m c + m s ){m 2 Da - [p + q) 2 ) Js s - [p + q) 2 

where p\ 2 show the spectral densities of the higher resonances and the continuum in 
hadronic representation. These spectral densities are approximated by evoking the quark- 
hadron duality assumption, 

pU")=P?,2 D (s)9(s-s ), (11) 
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where, pi^. D {s) = -IraU9 CD {s) are spectral densities in QCD side and s is continuum 
threshold in D s channel. 

The correlation function in QCD side, U® CD (s) is calculated by expanding the T product 
of the currents in ([5]) in terms of the DA's of the r\ meson with increasing twist in deep 
Euclidean region, where (p + q) 2 <C 0. After contracting out the c quark pair, we obtain 



n M (p,g) = ij d 4 xe^(r]\s^(l- l5 )S c (x)(l-^)s(0)\0), (12) 

where, S c (x) is the full propagator of c quark. 

The light cone expansion of the quark propagator in the external gluon field is made 
in 151 ] . The propagator receives contributions from higher Fock states proportional to 
the condensates of the operators qGq, qGGq and qqqq. In the present work, we neglect 
contributions with two gluons as well as four quark operators due to the fact that their 
contributions are small [16(. In this approximation, the S c (x) is given as: 



SJx) 



+ 



d A k 

i 



-ikx 



ft + m c 
k 2 — m 2 



Ws 



d A k 



-ikx 



du 



1 ljk + m c 



2 ( m l-k^ G ^ Ux)a " U 



m 2 — k : 



ux ll G tll '{ux) r ) v 



(13) 



where G^ v is the gluonic field strength tensor and g s is the strong coupling constant. We 
can rewrite the Eq. ( |T2l) as: 



n M (p,?) = - A IdV 



Tr 7/J (l - l5 )S c (x)(l - 75)^ (vlsF'slO) 



(14) 



where P is the full set of the Dirac matrices, P = (/, 75, 7a , 70,75, (J a p)- As it is clear 
from Eq. (114)) . to proceed to calculate the theoretical side of the correlation function, we 
need to know the matrix elements of the nonlocal operators between vacuum and 11 meson 

n 

states. Up to twist-4, the rj meson DA's are defined as [17[ : 



iv(p)\Tf^5q\0) = -ifnPp \ due lupx \<p v {u) + ^m 2 x z A(u) 



-7;fv m n— 



iupx \inJv) + —r- 2 - 2 

16 

2 ^ fl due- iupx B(u), 







2 JV v pxJo 



(15) 



(v(p)\q(xhn7bg s G a p(vx)q(0)\0) = f v m 2 



px 



Pa 9Pfi 



X/3PtM 

px 



x / Va^ ± (ai)e- ipx{ai+ua3) 
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(v(p)\Q(x)g a G a p(vx)'y lt g(Q)\0) = if v m: 

x / ' Vai<p±(ai)e 



(16) 

B\9an ~ Pal 9 fin 



px ) \ px 



+if ri m 2 ^-(p a x f3 - p p x a ) J Vai(p\\{ai)e 



2 EH(r, rr.„-r>„rr. \ I Tin, -mJ n> .\ p -ipx{<*i+ua 3 ) ^ 

(17) 



where, G^ u = \e^ vu \G aX and Voti = daida2da 3 5(l — a\ — a 2 — oi 3 ). Since we use the chiral 
current, the twist-3 wave functions do not give any contribution. In Eqs. (fT5l)-( TT7l) . the 
(Pr){u) is the leading twist-2, A(u) and part of B{u) are two particle twist-4, (p^ai), ip±(aii), 
<p\\{ati) and <£>_!_ ( a i) are three particle twist-4 DA's. Here we should stress that using the 
identity, 

l^af) = i{g^alf} ~ 9rfla) + Wp7 P 75 , (18) 

and due to the parity invariance of strong interactions, the matrix element, 

( V (p)\s^G a P(ux)o- a ps\0) = 0, (19) 

and has no contribution. For extracting the QCD or theoretical side of the correlation 
function, we insert the expression of the charm quark full propagator as well as the DA's 
of the r\ meson into Eq. (Till) and carry out the Fourier transformation. 

Now, we proceed to get the LCSR for our form factors equating the coefficients of the 
corresponding p^ and g M structures from both phenomenological and QCD sides of the 
correlation function and applying Borel transform with respect to the variable (p + q) 2 in 
order to suppress the contributions of the higher states and continuum as well as eliminate 
the subtraction terms. As a result, the following sum rules for the form factors f+ B ^ v and 

j-Ds^V + fDs^v are obtained: 

f u.->„/ 2 n z : m l m Vn ^ff f l du( 2^(m) 3A(u) _ m 2 c A(u) \ 

J+ W ) r> 2 £ 6 M 1 /. „. A..-\/T<> n„.0*/rA 6 M 



2 ™dJd s [Js u \ m 2 AuM 2 2u 2 M A 

f 1 , f u , B(u) 
+2 / du dt—^re~^ r 
Js Jo tM 2 



s J 1 k 2 M 2 



+4m 2 



du 



v / k dt ^M + n («0 - 2y±(a0 - 2^|K) e ^l (2Q) 
o t 2 M 4 1 



m 



W v fr, 



m2 Ds 

e m 2 



2 du dt-^-e^ 
J 6 Jo t 2 M 2 



(t) 



- S (t) 



o 



t 3 M 4 



-e «" 2 



(21) 



where, M 2 is the Borel parameter and, 

m 2 — g 2 x + mlxx 
s(x) = — — 

x = 1 — X 

k 

5 



«i + ua 3 , 
1 

' 'm 2 + g 2 - s ) + 



2m 2 



s — m 2 — g 2 ) 2 — 4m 2 (g 2 



(22) 



III. NUMERICAL ANALYSIS 



In this section, we numerical analyze the form factors, f± 



(g 2 ) and calculate 



branching fractions of the D s — > (j],T}')lv decays and their ratio. We also compare the 
results of the considered observables with predictions of the other nonperturbative ap- 
proaches as well as existing experimental data. As we mentioned before, using Eq. (j3J), the 
transition form factors of D s — > rj'lu decay are calculated by the help of the transition form 
factors of D s — > rjlu decay easily. Hence, we will discuss only the f± s ^ v {q 2 ) form factors. 
From the LCSR for these form factors, it follows that the main input parameters are the 
DA's of the rj meson. The explicit expressions of the wave functions, <f v {u), A(u), B(u] 
and (p\\(azi), (p±(aii), <pn(ai), and <p±(ai) as well as related parameters are given as [17] : 



<PvW 



6uu (l + a v 2 C£{2u- 1] 
120aia 2 tt3 (^oo + v w (Sa 3 - 1)) 
120aia 2 a 3 (cti (a 2 — «i)) , 



-30a 2 



h 0Q (l - a 3 ) + h 01 (a 3 (l - a 3 ) - 6a 2 «i) + h w (a 3 (l - a 3 ) 



/ 2 2\ \ 

-{a l + a 2 }) 



<P±(ati) = 30a 2 (ai - a 2 ) 



h 00 + h 01 a 3 + -h w (5a 3 - 3) 



B{u) 

9n( u ) 
A{u) 



g v (u) - (pr,{u), 

g cl(2u - 1) + g 2 cl(2u - 1) + 9A p\{2u - 1) 
6uu 



16 24 „ 20 / 1 1 7 10 , - 

T5 + 35°= + 20 " 3 + + I" 15 + 16 - 27*<" 3 " 27* 1 ° 2 (2 " " '< 



+ 



11 



210 
1? 



4 

135' 



V73W3 C 4 2 {2u - 1) 



+ ( _ y°2 + 2I774W4J [2m 3 (10 - 15w + Qu 2 ) \nu 
+ 2w 3 (10 - 15w + 6w 2 )lnw + ww(2 + 13mm) 

where C*(x) are the Gegenbauer polynomials, 

1 



(23) 



^-00 

fio 

^01 

hio 
9o 

92 
94 



21 



9 

^4^4 - —02, 



21 



7 

7 3 
4^4 + — a\, 



I^Wi - — ai, 



18 _ 20 
1 + y a 2 + 6O7/3 + y?74, 



9 

28 



- 6773^3. 



(24) 



The constants in the Eqs. ( 12"3"|) and (|24l) were calculated at the renormalization scale « = 
1 GeV 2 using QCD sum rules and are given as al = 0.2, i] 3 = 0.013, 7/4 = 0.5, W3 = —3 
and u>4 = 0.2. 

The values of the other input parameters appearing in sum rules for form factors are: 
qauark masses at the scale of about 1 GeV m s = 0.14 GeV, m c = 1.3 GeV 18|, meson 



masses m„ = 0.5478 GeV, = 0.9578 GeV, m Ds = 1.9685 GeV, V c , = 1.023 ± 0.036 



19] 



and f Ds = (0.274 ± 0.013 ± 0.007) GeV |20|. 

The sum rules for form factors also contain two auxiliary parameters, so and M 2 . The 
continuum threshold is not totally arbitrary but it depends on the energy of the first excited 
state. We choose, s = (6.5 ± 0.5) GeV 2 (see also [211]) . Now, we are looking for a working 
region for M 2 , where according to sum rules philosophy, our numerical results be stable 
for a given continuum threshold sq. The working region for the Borel mass parameter 



8 



is determined requiring that not only contributions of the higher states and continuum 
effectively suppress, but also contributions of the DA's with higher twists are small. Our 
numerical analysis shows that the suitable region is: 2.5 GeV 2 < M 2 < 3.5 GeV 2 . The 
dependence of the form factors f+ s ^ v and f Da ^ v on M 2 are shown in Fig. [TJ This figure 




FIG. 1: The dependence of the form factors on M 2 . The dashed, solid and dashed-dotted lines 
correspond to the so = 5.5, sq = 6 and sq = 6.5, respectively. 

shows that the form factors weakly depend on the Borel mass parameter in its working 
region. 

Now, we proceed to find the q 2 dependence of the form factors. It should be stressed 
that in the region, q 2 > 1.4 GeV 2 the applicability of the LCSR is problematic. In order to 
extend our results to the whole physical region, we look for a parametrization of the form 
factors such that in the region, < q 2 < 1.4 GeV 2 , the results obtained from the above- 
mentioned parametrization coincide well with the light cone QCD sum rules predictions. 
The most simple parametrization of the q 2 dependence of the form factors is expressed in 
terms of three parameters in the following form: 

/±(0) 



f±(q 2 



(25) 



1 — aq + /3q 2 

where, q = q 2 /m 2 Dg . The values of the parameters, /± S_H? (0), a and (3 are given in Table HI 
This Table also contains predictions of the light-front quark model (LFQM) for f+ e ~* v (0) 
for two sets (for details see |23(). The errors presented in this Table are due to variation of 
the continuum threshold So, variation of the Borel parameter M 2 , and uncertainties coming 
from the DA's and other input parameters. 

The dependence of the form factors, f+{q 2 ) and f-{q 2 ) for D s ^ rj on q 2 extracted from 
the fit parametrization are shown in Fig. (T5]). This figure also contains the form factors 
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TABLE I: Parameters appearing in the fit function for form factors of D s — > rj in two approaches. 



Model 


f D ^(0) 


a 


P 


This Work(LCSR) 


-0.44 ±0.13 


2.05 ±0.65 


1.08 ±0.35 




/+ s ^(0) 


a 


P 


This Work (LCSR) 


0.45 ±0.14 


1.96 ±0.63 


1.12 ±0.36 


LFQM(I)[23j 


0.50 


1.17 


0.34 


LFQM(II)[22] 


0.48 


1.11 


0.25 



obtained directly from our sum rules in reliable region. We see that, the aforementioned 
fit parametrization describe our form factors well. The values of the f+ s ^ <yr, ' n \q 2 ) form 





FIG. 2: The dependence of the form factors of D s — > r\ on q 2 . The circle points correspond to 
the values obtained directly from sum rules and the solid lines belong to the fit parametrization 
of the form factors. 

factors at q 2 = extracted from fit parametrization and using Eq. (j4j) are shown in Table 
HT1 Note that for massless leptons, the form factors, f^ s ^ v,r] \q 2 ) do not contribute to the 



decay rate formula, so we present only the /. 



(q 2 ) in this Table. For comparison, the 



predictions of the other approaches are also presented in this Table. From this Table, we 
see a good consistency among the results predicted by different approaches. 

Now, we would like to evaluate the branching ratios for the considered decays. Using 
the parametrization of the transition matrix elements in terms of form factors, in massless 
lepton case, we get: 

G 2 F \V CS \ 2 



1927r 3 m 3 



2\2 a 2 2 ' I t 

Q ) -4m D m M) \f_ 



3 / 2 |^->^ (g 2 )|2 



(26) 
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TABLE II: The f® a ^( ri ' ri \q 2 ) form factors at q 2 = in different approaches: this work 
(LCSR), three-point QCD sum rules (3PSR) and LFQM. Our results for /+ S_5>T? correspond to 
^ = 39.7°(41.5°). 



Form factor 


This work (LCSR) 


3PSR[22j 


LFQM(I)[23] 


LFQM (II) [23] 


/£•->) 


0.45 ±0.14 


0.50 ± 0.04 


0.50 


0.48 




0.55 ±0.18(0.51 ±0.16) 




0.62 


0.60 



where Gp is the Fermi constant. Integrating Eq. f l26|) over q 2 in the whole physical re- 
gion and using the total mean lifetime, td s = (0.5 ± 0.007) ps [l9|, the branching ratios 
of the D s — > (r},rf)lv decays are obtained as presented in Table fTTTl This Table also in- 



TABLE III: The branching ratios in different models and experiment. Our values correspond to 
39.7°(41.5°). 



Mode 


This work 


3PSP[22] LFQM(I)[23j 


LFQM(II)[23J 


EXP [19] 


Br(L> s 


— > rjlu) 


x 10 2 


3.15 ±0.97 


2.3 ±0.4 2.42 


2.25 


2.9 ±0.6 


Br(D s 


-»• rflv) 


x 10 2 


0.97 ±0.38(0.84 ±0.34) 


1.0 ±0.2 0.95 


0.91 


1.02 ±0.33 



eludes a comparison of our results and predictions of the other nonperturbative approaches 
including the LFQM and 3PSR and experimental values [l9|. From this Table, we see a 
good consistency between our results and predictions of the different approaches especially 
experimental data. 

At the end of this section, we would like to compare also the ratio: Rn = Br(D s -»r? iu) 
in Table [IV] for different approaches as well as experimental value. This Table also depicts 
a good consistency among the values, specially between our prediction with if = 39.7° and 
experimental value. This can be considered as a good test for correctness of the considered 
internal structure for the D s meson as well as the mixing angle between rj and rj states. 
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TABLE IV: The Rd s with respect to mixing angle, <p for different models and experimental value. 



Model Angle (ip°) R Ds 

This work (LCSR) 39.7°(41.5°) 0.32 ± 0.02(0.27 ± 0.01) 

3PSR[22J 40° 0.44 ± 0.01 

LFQM(I)[23J 39° 0.39 

LFQM(II)[23J 39° 0.41 

EXP [12] - 0.35 ±0.12 
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